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ABSTRACT 

We examine Hamiltonian formalism on Euclidean Snyder space. The latter corresponds 
to a lattice in the quantum theory. For any given dynamical system, it may not be possible 
to identify time with a real number parametrizing the evolution in the quantum theory. The 
alternative requires the introduction of a dynamical time operator. We obtain the dynamical 
time operator for the relativistic (nonrelativistic) particle, and use it to construct the generators 
of Poincare (Galilei) group on Snyder space. 
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1 Introduction 



Snyder showed long ago that a spatial lattice can be consistent with the continuous symmetries 
of space-time through the construction of a covariant noncommutative algebra. [1] The algebra 
is generated by two Lorentz vectors: the time-space coordinate x M , and the energy- momentum 
p^, jj, = 0, 1, 2, 3. The spatial lattice of Snyder, or 'Snyder space', results from the fact that the 
spatial coordinates x j, i = 1,2, 3, have discrete spectra. It differs from a classical lattice because 
only one coordinate can be determined in a measurement due to the noncommutativity. In a 
previous article, [2] we found that two distinct representations of the Snyder algebra are possible. 
They are characterized by an SU(2) quantum number, which takes on integer values for one 
representation and half-integer values for the other. Continuous symmetry transformations are 
unitarily implementable on the lattice. The continuous transformations examined in [2] were 
generated by the momenta and angular momenta, and are associated with the translation* 
and rotation group, respectively. In this article we shall demonstrate, among other things, 
how to write down the Hermitean generators for the full Poincare group (Galilei group). This 
is facilitated by examining the Hamiltonian dynamics of a relativistic particle (nonrelativistic 
particle) on Snyder space. 

Covariance is problematic in the standard Hamiltonian description of a relativistic particle 
since one treats time differently from the spatial coordinates. It is then not surprising that 
the Lorentz covariant algebra of Snyder is inconsistent with the Hamiltonian description of a 
relativistic particle. Various tactics can be taken to get around this issue, and they have been 
employed in previous articles. These include introducing additional degrees of freedom [3], such 
as an additional time variable [4], or modifying the particle action [5]. The approach adopted in 
this article does not require either introducing additional degrees of freedom or deforming the 
dynamics. Instead, as in the standard Hamiltonian description, we can drop the insistence on 
covariance, and still have a consistent relativistic particle dynamics. The latter is written down 
on the Euclidean subalgebra of the Snyder algebra, generated by the spatial components of the 
four vectors, i.e., Xi and pi.§ Although, p° and x°, no longer appear as independent generators 
of the algebra in this approach, the notion of energy and time re-appear with the introduction 
of particle dynamics. One can then derive algebraic properties for these quantities. This is 
straightforward for the operator p° associated with the energy, once we define it by the mass 
shell condition. 

Different approaches can be taken with regard to the time. If one insists, as usual, that time 
is a real parameter associated with the evolution of the system, then p° cannot be the particle 
Hamiltonian on Snyder space because it will not in general generate time evolution. Rather, 
the Hamiltonian is some function of p°. This then implies nonstandard energy- momentum 
dispersion relations, a feature in common with double special relativity [7]. On the other hand, 

was shown in [2], they are not translations from one point on the lattice to another, but rather are 
translations in a particular continuous basis of the Hilbert space. 

§ Some aspects of nonrelativistic dynamics on this space were examined previously in [6]. 
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as we argue here, nonstandard energy-momentum dispersion relations need not be a signature 
of Snyder space. One can retain the conventional energy-momentum dispersion relation upon 
re-introducing a time operator. The time operator in this case is not an independent generator 
of the algebra, but rather a function of the operators Xj and pi, along with the evolution 
parameter denoted by A. For this reason we refer to it as a dynamical time operator. The 
operator is only defined up to constants of the motion, which are not in general central in the 
algebra. For several examples we determine the definition of the time by requiring it to be 
a gauge fixing condition in a reparametrization invariant formulation of the dynamics. The 
choice for the gauge fixing is such that the Euclidean Snyder algebra is realized by the Dirac 
brackets for the constrained Hamiltonian system. The examples considered in this article 
are the relativistic and nonrelativistic free particle, along with one-dimensional conservative 
dynamics. The relativistic (nonrelativistic) free particle action that we start from is Lorentz 
(Galilei) invariant, as well as reparametrization invariant. We are then able to construct 
Lorentz (Galilei) boost operators from the Euclidean Snyder algebra which satisfy the Poincare 
(Galilei) algebra. 

The outline of this article is as follows. In Section 2 we show that the covariant Snyder 
algebra is inconsistent with the mass shell condition. We then argue that dynamics can be 
consistently formulated on the Euclidean subalgebra. In Section 3, we write down the symplec- 
tic two-form associated with the Euclidean Snyder algebra and develop Hamiltonian dynamics 
on the resulting phase space. Time is treated as a commuting number in this section, while 
dynamical time operators are constructed in Section 4. We do this for the three examples of 
the nonrelativistic and free relativistic particle, and the one-dimensional conservative system. 
Momentum-dependent eigenfunctions of the time operator, along with space-time symmetry 
generators, are obtained for the free particle examples. Concluding remarks are made in Section 
5. 



2 Covariance lost 

The Lorentz covariant algebra of Snyder is defined by the commutation relations [1] 



^2 (^IJ-Pv x vPn) 

Ppi>v 



[x^Pu\ = l[V^u+ A2 

\p»,Pu] = (2.1) 

where fx, v, ... = 0, 1, 2, 3 and we choose [t/ui/]= diag(— 1, 1, 1, 1) for the Minkowski metric. A / 
is the deformation parameter, which has units of energy. The canonical commutation relations 
of Minkowski space-time coordinates with the four- momenta are recovered in the limit A — > oo. 
The algebra (2.1) was originally obtained starting from the de Sitter group acting on de Sitter 
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momentum space, and then projecting to Minkowski space-time. The projection involves a 
particular coordinatization of de Sitter space by the four-momentum p„, and an identification 
of four of the de Sitter group generators with the space-time coordinates x^. Following that 
approach one arrives at an upper bound on the mass, —p^p^ < A 2 , but this bound is not 
required if we instead simply postulate the commutation relations (2.1) from the start. 

The canonical commutation relations of Minkowski space-time coordinates and the four- 
momenta (corresponding to the limit A — > oo) are inconsistent with the mass shell condition 
P^Pn + m 2 = defining a relativistic free particle. This is since p^p^ is not central in the 
algebra. It is then not surprising that the algebraic relations (2.1) are also inconsistent with 
the mass shell condition. This is since 

[W"A,] = 2^(1 + ^) (2.2) 

Just as one should not insist on Lorentz covariance when writing down the Hamiltonian for- 
malism for a relativistic particle using the canonical commutation relations, one should not 
insist on Lorentz covariance when writing down relativistic particle dynamics on Snyder space. 

Instead, Hamiltonian dynamics for a particle moving in three dimensions is consistently 
written down on a six-dimensional phase space. By Darboux's theorem we can take the latter 
to be the Euclidean subalgebra of the Snyder algebra, at least classically. At the quantum level, 
the algebra is generated by the three spatial components of the coordinates and momenta, X{ 
and pi, i = 1, 2, 3, and defined by the commutation relations 

Pipj 



[xi,pj\ = i\6ij+ A2 

\pi,Pj] = (2.3) 

where Lj = eijkXjpk are the angular momenta. The three-momenta pi are simultaneously diago- 
nalizable and have continuous eigenvalues pi . Snyder gave a differential operator representation 
for x-i acting on the space of complex functions {(j), ip, ...} of the momentum p = {pi-,P2-,Pz)- It 
is:1 

^.Mi (2 4) 

These operators are symmetric for the scalar product 

(0,V)= / dn{p) m*1>(p) , (2-5) 



using the measure 

dn(p) = d \ (2.6) 
(1 + g )2 



^This was generalized to a one parameter family of representations in [2]. 



4 



Whereas Snyder's Lorentz covariant algebra (2.1) possesses the independent operators x° 
and p°, associated respectively with the time and energy, from (2.3) we can define analogous 
operators as functions of x % and p l using particle dynamics. For the relativistic free particle, 
the energy is naturally defined by the mass shell condition 

P° = VPiPi + m2 (2-7) 

Then 

fc.j?] + [p l ,p°}=0, (2.8) 

which differs from the commutation relations for p° postulated in (2.1). In the limit A — > oo, 
(2.8) implies the usual Heisenberg equation of motion for a relativistic particle upon identifying 
p° with the Hamiltonian. 

As mentioned in the introduction, the notion of time for a relativistic particle can be 
addressed in different ways. It can be regarded as the evolution parameter associated with the 
relevant Hamiltonian H t , or it can be an operator defined in terms of some evolution generated 
by the zeroth-component of the four- momentum p°. Both of these approaches coincide, when 
one uses the Heisenberg commutation relations, i.e., A — > oo, since then Ht = p°. However, this 
is not true in general, and specifically for commutation relations (2.3). In the former approach, 
we avoid dealing with a noncommuting time operator. Then Ht is some deformation of p , 
which we examine in the following section. The latter approach re-introduces a time operator, 
which here is a function of the fundamental operators and pi, as well as the evolution 
parameter for the system, which we denote by A. As stated in the introduction, the function 
may be interpreted as a gauge fixing condition starting from a reparametrization invariant 
action. The commutation relations for this dynamically defined time operator can then be 
derived using (2.3). The advantage of this approach is that one can retain j5° as the evolution 
operator for the relativistic particle [or fP/(2m) as the Hamiltonian for the nonrelativistic 
particle] and therefore preserve the standard energy-momentum dispersion relations. This 
appraoch is examined in section four. 



3 Hamiltonian formalism on Snyder space 
3.1 General formalism 

Before discussing quantum dynamics on the Snyder algebra, we examine classical Hamiltonian 
dynamics on the associated classical phase space. The phase space is spanned by X\ and pi, 
corresponding to the classical analogues of X{ and pi. The symplectic two-form for the system 
is given by 

co = dxi A dpi - -d(xiPi) A d\n(j? + A 2 ) (3.1) 

To see this, we can invert the symplectic matrix to obtain the associated Poisson brackets. 
The inverse exists because the determinant is (1 +pP/A 2 )~ 2 , which is nonvanishing for finite 
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momentum. The resulting Poisson brackets are 

1 

A 2 " 



{Xi,Xj} — 2 £ijkL>k > L{ — £ijk%jPk 



{ x i}Pj} — $ij + 

{Pi,Pj} = 0, (3.2) 

which are the classical analogues of the commutation relations (2.3). The second term in (3.1) 
drops out in the limit A — > oo, and we recover the canonical symplectic two-form 

ojq = dx,i A dpi (3-3) 

More generally, the Darboux map from (3.1) to the canonical symplectic two-form (here written 
as dqi A dpi) is given by 

Xi — > qi = Xi — ^ _^p2 ' 

while the momenta pi are invariant under the map. Thus the variables qi are canonically 
conjugate to p{. 

For dynamics on the phase space spanned by X{ and pi, we introduce a Hamiltonian H. 
It generates evolution in some parameter A. In general, A need not be the time t, but rather 
some monotonically increasing function of the time. The Hamilton equations of motion 

i& x oj = dH (3.5) 

relate H to the dynamical vector field 

A = dxj__d_ + djH_d_ , 3g s 
A dX dxi dX dpi 

From the symplectic two-form (3.1), one gets the following equations for the first derivatives 
of H, 

dH d Pl , ^(f) Pi 



+ 



dxi dX 2(p 2 + A 2 ) 

dH dX, &(j?) Xi -2&(x-p) Pi 

d Pi dX 2(p2 + A 2 ) 1 ' ' 

Consistency relations can be obtained by computing second derivatives of H . 

Aspects of Hamiltonian dynamics on the above phase space were examined previously in 
[6], [10], along with the equations of motion following from sample Hamiltonians, in particular, 
the harmonic oscillator Hamiltonian 

Hho = ^~ + \mu?x 2 (3.8) 
2m 2 
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The evolution parameter A was equal to the time in these works. With such an identification, 
the dynamics generated by H on the phase space associated with symplectic two-form (3.1) is, 
in general, inequivalent to the dynamics generated by the same Hamiltonian on canonical phase 
space, the latter defined by the symplectic two-form (3.3). Then, for example, along with 
(3.1) does not describe a simple harmonic oscillator, but rather a deformed oscillator. Such a 
system was examined in [6]. In this article our aim is not to deform, but rather, to preserve 
the dynamics of any given system while passing from the canonical symplectic two-form (3.3) 
to (3.1). 

To identify the equations of motion (3.5) with some given system, we must specify both 
H and A. Say that the dynamics is Hamiltonian with respect to the canonical symplectic 
two-form (3.3). That is, it is defined by the symplectic two- form (3.3) and some Hamiltonian 
Hq, with the latter generating evolution in the time t.H Thus the dynamics is given by the 
standard Hamilton equations 

iA t ^o = dH , (3.9) 
where A t = |^A A = ^f^f- + ^ff g§r, or in terms of components, 

8Hq dpi 8Hq dxi 



dxi dt dpi dt 

Substituting (3.10) into (3.7) gives 

d\ dH 8Hq piPj dHo 

dt dxi dxi p 1 + A 2 dxj 

d\dH 8H If, ,9H dH 



(3.10) 



dt dpi dpi p 1 + A 2 



( , ,9H dH \ 

\J XiPj - XjPi )— +PiPj —) (3-H) 



The conditions state that H generates the same dynamics on the phase space manifold with 
symplectic two-form (3.1), as Ho generates on the phase space manifold with the canonical 
symplectic two-form (3.3). The equations of motion (3.11) are valid in any number of dimen- 
sions. For a given Hq, one can search for consistent solutions to H and 4| as functions of the 
phase space variables x% and pi, and also A. These solutions, if they exist, may not be unique. 
Below we shall obtain solutions in two simple cases: the case of one-dimensional motion and 
the case of free particle motion. 



3.2 One-dimensional systems 

We first consider the case of motion in one dimension, by which we mean a two-dimensional 
phase space spanned by x and p. Then Eqs. (3.11) simplify to 

dX dH A 2 dH 

dt dx p 2 + A 2 dx 



"More generally, it may be possible to have is a solution to (3.7) even when no Hamiltonian description with 
respect to the canonical symplectic two-form u>o is possible. 
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d\8H_ _ A 2 8H 

dt dp p 2 + a 2 dp { ' 

They are satisfied upon writing 

H = f(H ) ^ = (l + j£) f'(H ) , (3.13) 

where /{Hq) is a function of Hq, and f'{Ho) is its derivative. It remains to specify the evolution 
parameter A. Here we should require that ^ is everywhere positive for all physical particle 
trajectories. Of course, this is not a concern if we simply identify A with the time t. In that 
case, from (3.13), we would need f'(Ho) = (1 + p 2 /A 2 )~ 1 . This can be satisfied when Hq 
depends only on p, but it is not in general valid for an arbitrary x-dependent Hamiltonian, 
such as the case with a harmonic oscillator. 

Alternatively, instead of trying to identify t with A, one can set H = H , and then solve 
for A as a function of x, p and t. This is always possible on the two-dimensional phase space. 
In that case, one has 

dt v 2 

H - H " 5A = 1 + 7P- < 3 ' 14 > 

and it follows that ^ is everywhere positive.** To obtain the relation between t and A we 
should integrate (3.14) on the space of solutions: x = x^f 2 ' '"(i), p = plof 2, "'(t), where £i,£2, ••• 
label the solutions and correspond to constants of motion for the system. Then 

A = A 2 / f- , (3.15) 

J A 2 + \pt^-(t>)] 2 

from which one can in principle solve for t as a function of A and the phase space variables. Thus 
t is dynamically determined, and in general will have nontrivial Poisson brackets. However, 
the latter are not uniquely determined, since one can add a function of the constants of the 
motion £1,^2, ••• to (3.15). The constants of motion are not, in general, central in the algebra, 
so the commutation relations for the time operator depends on this function. In section 4, we 
fix this arbitrariness by demanding that the symplectic structure (3.1) be derived starting from 
a particle action. We require the action to be reparametrization invariant. The Snyder algebra 
is then realized by the resulting Dirac brackets for the system upon fixing a gauge condition. 
The gauge fixing condition then defines the time in terms of the evolution parameter. (This 
procedure was previously carried out for the case of the relativistic free particle in [8] . Here 
we also apply it to the nonrelativistic particle and one-dimensional conservative systems. ) 



3.3 Free particles 

The problem of finding solutions to (3.11) in more than one dimension is more complicated - 
with the exception of the free particle case, which we now discuss. Here we assume that Hq is 



**We assume that A 2 > 0. 
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independent of Xj. (3.11) then implies that H also only depends on the momenta (and possibly 
A) and the dependence is such that 

— —-fx PiPj \ dH ° (*-iR\ 

dt d Pi ~ V j p 2 + a 2 ) dp, [6 bj 

Upon insisting on rotational invariance, i.e., Hq is a function only of p 2 , and the equations 
simplify further to 

dXOH A 2 8H (; , i?) 



dt dpi p 2 + A 2 dpi 

The solutions are once again (3.13). It again remains to fix H and A. Here given Hq one can 
consistently set t = A and solve for H as a function of p 2 . We do this explicitly below for the 
nonrelativistic and relativistic free particle. 

Starting with the nonrelativistic free particle, with Hq = we get 

A 2 / f\ 

t = X H = H t = -ln(l + ^) (3.18) 

up to constants of integration. Regarding Hi as the energy, (3.18) implies a non-standard 
energy-momentum dispersion relation. Upon expanding Hi in 1/A, the standard nonrelativistic 
kinetic energy appears at zeroth order and the leading correction is fourth order in momentum, 

H t = ^—- ^L* + • • • (3.19) 
2m 4mA 2 v ; 

For the case of the relativistic free particle, we can choose H = p°, where 

p° = ^p^ + m z (3.20) 

It leads to the equations of motion dx/dt = p/p°. Here t should be interpreted as the zero 
component of the space-time four vector. Up to constants of integration, the Hamiltonian H 
which generates the motion on the phase space associated with symplectic two-form (3.1) is 

( /It i ta 11 " 1 rrr i » m < A 

VA 2 -rrt 2 VA 2 -m 2 

t = A H = H t = I (3.21) 

- / k i a 2 tanlT 1 , \ 2 , m > A 

Vm -A 2 Vm -A 2 

(This result appeared previously in [8].) Regarding H t as the energy, we once again get a 
nonstandard dispersion relation for the energy-momentum for finite A. One recovers p° from 
(3.21) in the limit A — > oo. The leading 1/A 2 correction is 



rn 



2 



(2m 2 -p 2 ), (3.22) 



3A 2 

which includes a |^ correction to the rest mass energy. 

For both of the above free particle solutions, time remains a commutative parameter in 
the quantum theory. It, along with operator analogue H t of the Hamiltonian H t , defines the 
quantum evolution operator 

U(t, 0) = e~ iAtt (3.23) 
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From the evolution operator one recovers the usual equations for the expectation value of the 
operators. For the nonrelativistic particle, one gets 

±<Xi>=-i< [xi, H t ] >= ± < m >= , (3.24) 

at m at 

while for the relativistic particle 

j t <x i >=-i<[x i ,H t }>=<^> ^<pi>=0, (3.25) 

where j5° is defined in (2.7). 

Unlike in the above, generic solutions to (3.11) will have ^ / 1, and so then time cannot 
be identified with the evolution parameter, either in the classical or the quantum theory. In 
that case, time is a function of the position and momenta, as well as A, in the classical theory, 
while it is an operator in the corresponding quantum theory. We examine this possibility, 
specifically for the case H = Hq, in the next section. 



4 Dynamical time operators 

In this section we choose H to be the Hamiltonian Hq which generates evolution on the canon- 
ical phase space. By so doing, we avoid nonstandard energy-momentum dispersion relations, 
such as in (3.18) and (3.21). We obtain the time t as a function of x, p and A in the clas- 
sical theory for three examples, the free nonrelativistic particle, the free relativistic particle 
and one-dimensional conservative systems. We do this by demanding that it corresponds to a 
gauge fixing condition for a reparametrization invariant action. 



4.1 Nonrelativistic free particle 

Since the momentum p is constant for free particle trajectories, (3.14) implies that t is a 
rescaling of A by a constant factor. From (3.15) one gets 

along any trajectory. C is some function of the constants of motion, which here include the 
Galilei boosts 

G = tp-mx, (4.2) 
as well as p and L. Below we shall show that if we make the choice C = p- G, and hence 

mx ■ p ... 

that we recover the full Galilei symmetry for the nonrelativistic free particle and, in particular, 
Gi satisfy the algebra of Galilei boosts. 
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Galilei invariance is, of course, insured if one starts with the usual free particle action 
^nr = f dt ^(^^j ■ The action can be rewritten so that it possesses a reparametrization 
symmetry. Following [11], one can preserve the Galilei invariance, while introducing a gauge 
degree of freedom i(A) associated with the time [where A is an evolution parameter with no a 
priori relation to t]. The corresponding action is 



/ 



f I , (4.4) 



the dot denoting differentiation in A. If pi denotes the momenta conjugate to %i and irt denotes 
the momentum conjugate to t one gets 

•_,2 

Xi m x 

Pi = mj 7r t = -— -p-, (4.5) 

and hence the Hamiltonian constraint 

xjr 1 = f + 2mvr t ps , (4.6) 

generating the reparametrization symmetry A — > A' (A). The symmetry is eliminated upon 
identifying A with the time t, in which case (4.4) reduces to 5„ r . Here let us instead use (4.3) 
to relate A to the time t. That is, we impose the gauge constraint 

¥ 2 = t -A-^«0 (4.7) 

It is straightforward to show that the resulting Dirac brackets for Xi and pi are (3.2), yielding 
the Euclidean Snyder algebra. The Dirac brackets of the time t with Xi and pi are 

{xi, t} DB = -p- \5ij - A2 + j {pi,t} D B = — (4.8) 

Since the action (4.4) is Galilei invariant and leads to the Snyder algebra when imposing the 
gauge fixing (4.7), Galilei invariance must follow when defining the dynamical time for the 
Snyder algebra according to (4.3). 

Next we promote the classical time t to a Hermitean operator £ nr . Up to operator ordering 
ambiguities, it is 

1/ -2 -2 \ 

tnr = A + - [XiPif nT (p ) + /nr(p )piXij , (4.9) 

where 

Ur(p) = m(A 2 +f)- 1 (4.10) 

Since the time operator depends on Xi and pi, it is dynamically determined. From (2.3), we 
derive the following commutation properties for the dynamical time operator: 

if ->2 -*2 \ 

[Xi,t m ] = ^2 [mXi - XjPjPif nr (p ) - /nr(p )piPjXjj 
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IfJl 

[Pi,tnr] = -^Pi, (4.11) 

which are the quantum analogues of (4.8). It follows that t nr is rotationally invariant, [Lj, t nT ] = 
0. The commutation relations (4.11) differ from those postulated by Snyder for x° in (2.1), 
which is obvious since Lorentz covariance is not present in (4.11). 

The quantum analogue of (4.2) is 

Gi = ^(i m ft + Pii nr ) - mxi (4.12) 

Gi satisfy the algebra of Galilei boosts 

[Gi,pj] = -imSij 

[Gi,Lj] = ieijkGk 

[Gi,H ] = -ifr (4.13) 

Moreover, the full Galilei algebra is recovered upon including the generators pi, Li and the 
Hamiltonian 

h = M (4 . 14) 

Here, in addition to (4.13), one has the commutation relations 

[Li,Pj] = ieijkPk 

[Li,Lj] = ie ijk L k , (4.15) 

and \pi,Pj] = [pi, Ho] = [Li, Ho] = [Gi,Gj] = 0. The Galilei group acts unitarily in the 
quantum theory. The unitary actions of pi and Li were given explicitly in [2], and using them 
one can construct the remaining Galilei transformations. This follows because the remaining 
generators Gi and Hq are written in terms of Xi and pi which have a well-defined action on 
the states. 

As i nr is noncommuting, care must be taken in writing down a time-evolution operator in 
the quantum system. It is constructed using the Hamiltonian (4.14) according to 

U(\,0) = e~ iAoX (4.16) 

The evolution of the mean values of the space and time coordinates is determined from 

-a 



d r- tt i /Pi 

-2 



< o > - -/ < [:>-,. /:/ ] > - — (• + 



A. < i nr > = -i < [£ nr , H ] > + (^) = (l + ^ ) , (4-17) 
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and so the mean velocity is 

d < Xi> _ (m ( 1 + I 2 ") 



d < I in > 1 1 _|_ t 



(4.18) 

7?) 

in contrast to the usual result (3.24). The expression (4.18) reveals the unusual feature that the 
mean velocity is not linear in the momentum, and that the size and shape of the wavepacket 
can effect its mean velocity when one approaches the noncommutative scale A. To illustrate 
this point we can consider a Gaussian distribution in momentum space centered about p = 
(0,0,<p>), 

^(p)Ia=o = ^37^572 ( 1 + A2J eXp \ 2^ / (419) 

The factor 1 + was included to cancel out the momentum-dependent factor in the measure 
(2.6), and we assume a > 0. The width in momentum space Api is in any direction, 

i = 1,2,3. Using the differential representation for x\ in (2.4), we get that the wavepacket is 
centered about the origin in position space when A = 0, < Xi > | a=o = 0- For the width of the 
wavepacket in the £3 direction, we find 

(^^{(i + ^) 2 + ^( 1 + t1|^) + ^} , ,4,0, 

which reduces to the usual value of 1/ \^2a in the A — > 00 limit. The effect of noncommutativity 
is to increase the spread of the wavepacket in position space. This result is consistent with the 
modified Heisenberg uncertainty relation. From (4.20) one can check that 

Ax 3 A P3 > + > (4-21) 

which follows from the commutation relations (2.3), is satisfied at A = 0.^ At an arbitrary A, 
we get from (4.17) that 



ft We have assumed h = 1 throughout this article. For more discussion of such uncertainty relations, see 

[9], [10]. Here we get an additional time-energy uncertainty relation 

~ n 1 " 

Atnr AH > < H > , 

A 2 

■ -*2 

which follows from [t nr , Ho] = P ■ For the wavepacket (4.19), one gets < i nr > |x=o = and 



(AM^o) 2 = ^(3 + 



< p > 2 



T 2 

after applying the differential representation (2.4). This along with 



1 / 2 3(7 \ . ' <J 3(7 2 „ ,, 

<Ho>=^ (<P> 2 +— ), AH ° = ^— + 2 <P> 2 



is consistent with the inequality. 
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->2 



< t nr > = ( 1 + 



2 i 3^2 



^-(l + iH^tSf-)* (4.22) 

The resulting mean speed of the Gaussian wavepacket as a function of its width and of < p > 
is 

Q < x 3 > = Z + <p>^+A^ < g > 
d<i QT > "" 2 + ^|^ m 

For <j<<Aor|<p>|, the mean speed differs from <p>/m by a the factor of 1 + 
a 2 /(<p> 2 +A 2 ). 

Unlike the position operators, the time operator (4.9) has continuous eigenvalues. We 
denote the latter by r. Then using the differential representation (2.4), we can solve the 
eigenvalue equation 

tnr<Mp,A) = ri(p,A) (4.24) 
for the eigenfunctions in momentum space 4> T (p, A), evaluated at some A. We get 

Their normalization is such that 

dnip) <P T >(P, A)V T (p, A) = <5(r - r') (4.26) 



To see this, substitute (4.25) in the left-hand side and apply the integration measure (2.6) to 
get 

Then after making the change in integration variables from p to w = — lnp, this becomes 
^ dw e^ T '- T > = 5(t - t'). 



4.2 Relativistic free particle 

We now repeat the previous analysis for the relativistic case. A is again given by (4.1) along 
any trajectory, where C is some function of the constants of motion. From the equation of 
motion dx/dt = p/p°, the constants of motion are now p, L and 

K = tp — xp° (4.28) 

We again define p° using the mass shell condition (3.20). The full Poincare symmetry of the 
relativistic particle is recovered for the choice C = p- K, and hence 

p° x • p , 
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with the constants of motion generating the symmetry. Eq. (4.29), along with the equations 
of motion, leads to 

dx % pi ( p 2 \ dt p 2 / , x 

The second equation indicates that A is a momentum-dependent rescaling of the time t on any 
particle world line, as noted previously. Eq. (4.29) reduces to (4.3) in the nonrelativistic limit. 

Analogous to what happens in the nonrelativistic case, Poincare symmetry follows because 
one can write down a reparametrization invariant particle action, which yields the Snyder 
algebra when imposing (4.29) as a gauge fixing. [8] The action can be taken to be the standard 
one for a relativistic free particle action 

S rcl = -mJd\ yZ-x^xf , (4.31) 

with the dot again denoting differentiation in evolution parameter A [which we assume has no 
a priori relation to x ]. As is well known, the action leads to the mass shell constraint 

*i = P% + m 2 « , (4.32) 

in the Hamiltonian formalism, with p^ canonically conjugate to x^. The constraint generates 
the reparametrization symmetry, which may be eliminated by imposing an additional constraint 
on the time t = x°. It was shown in [8], that the Euclidean Snyder algebra (3.2), is realized 
by the Dirac brackets if one instead chooses the gauge fixing condition to be 

*2 = *-A-f^-|«0 (4.33) 
A z + p z 

Since the action (4.31) is Poincare invariant and leads to the Snyder algebra when imposing 
this gauge fixing (4.33), Poincare invariance must follow when defining the dynamical time for 
the Snyder algebra according to (4.29). Below we shall show explicitly that this symmetry is 
generated by p^, Li and Ki. 

Up to operator ordering ambiguities, the Hermitean operator i re \ associated with t in (4.29) 

is 

1 / -2 ->2 \ 

£ re l = A + - [XiPif rc l(p ) + / re i(p )piXij , (4.34) 

_2 

where the function f re i((p ) is given by 

/rel(| 2 ) = (A 2 +/ry, (4.35) 

and reduces to (4.10) in the nonrelativistic limit. t re \ is dynamically determined since it is a 
function of Xi and pi. It is inequivalent to the time operator x° appearing in (2.1). Unlike 
(2.1), the algebra generated by t re \, Xi,pi and p° [given in (2.7))] is not Lorentz covariant. From 
(2.3), we get the commutation relations 

\x- i ,1 - —(fPx- + x rP + A 2 ~ (P ) 2 ~ m 2 p^ pjfo A 2 -(p°) 2 -m 2 A 

ZA V A 2 + p P P A 2 + p / 
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[pVrel] = -J^PiPi (4-36) 

The second equation agrees with Snyder's commutator [pi,x ], while the remaining two equa- 
tions differ. The third equation follows from the second using (2.7). The time operator is once 
again rotationally invariant, [Li,t TC \} = 0. 

We recover the Poincare algebra, despite the lack of covariance. The quantum analogue of 
(4.28) is 

K i = 2 (*relPi + P»*rel ~ %iP° ~ P° x i) ( 4 -37) 
From it we obtain the standard commutation relations for the Lorentz boosts, 



[Ki,?] 


= -ipi 


[ki,Pj} 




[Ki,Lj] 




[Ki,kj] 


= —itijkL, 



(4.38) 

Upon including pi,p° and Lj, we obtain (4.15) and \pi,fij] = \pi,p°] = [Li,p°] = 0, and hence 
the Poincare algebra. The Hilbert space for the theory therefore carries unitary representations 
of Poincare group. 

Poincare transformations of space-time are implemented with the action of pi, po, Lj and 
Ki. To see this we need to insure that the gauge condition (4.29) is preserved while also 
performing a reparametrization A — > A' (A) in the action (4.31). One implements the Poincare 
transformations on space-time, as usual, using commutators with the generators, only here 
the position operators have discrete spectra. Below we simplify to the classical theory and 
show that infinitesimal space translations, time translations, rotations and Lorentz boosts are 
obtained by taking Poisson brackets, respectively, with the Poincare generators Pi, po, Li and 
Ki. 

1. Infinitesimal spatial translations leave the four- momenta invariant, while the space-time 
coordinates of the particle x % {\) and t(A) undergo variations 

6x l {\) = € l + ^-5X 6t(\) = -j-6\, (4.39) 

€i being infinitesimal parameters. It remains to find how 8 A depends on e l . Using (4.30), 

-2 -2 

8x l {A) = e i + ^{l + ^8A St{X)=(l + ^SX (4.40) 
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The infinitesimal reparametrization 5 A is constrained by the requirement that the gauge 
condition (4.29) is preserved by the variations (4.40), i.e., 



p° + Sx(\fj -p 

i -r ol{a ) = a -r 

This gives 



t + ft(A) = A+ \i +pi ' (4.41) 



It then follows that the infinitesimal variations (4.40) associated with spatial translations 
are obtained from Poisson brackets with pi, 



fa*(A) = (5 ij + P -^-)e j = {x\e-p} 



5t{\) = = {t,e-p} (4.43) 



5x*(A) = + J)<5A 5t(X) =e + [l + ^)s\, (4.44) 



2. Infinitesimal time translations also leave the four-momenta invariant, while 

where e° is infinitesimal and we again used (4.30). In this case the gauge condition (4.29) 
is preserved for 

SX = -e° , (4.45) 

and the infinitesimal variations (4.44) associated with time translations are obtained from 
Poisson brackets with po, 

sax) = -%(i+£y = 

5t(X) = -^e° = {t,e°p } (4.46) 

3. Infinitesimal rotations leave po and t invariant, while 

Sx l = e ijk rfx k 5pi = e ijk V J Pk , (4.47) 

where rj k are infinitesimal. The transformation involves no reparametrizations since the 
gauge condition (4.29) is rotationally invariant. As usual, (4.47) can be expressed in 
terms of Poisson brackets with , 

Sx i = {x\ff- L} 5 Pl = { Pi ,ff-L} (4.48) 



4. Finally, infinitesimal Lorentz boosts are of the form 

„^ 1 + 



-Q. -2 

5x l (X) = co°H(X) + Pi(l + ^jSX 6t(X) = w°V(A) + (l + ^)SX 
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S Pi = -lu^po 5 Po = -u» l Pi , (4.49) 

where oj 1 ° are infinitesimal. In order for the gauge condition (4.29) to be preserved, one 
now needs 

SX = ^^{(A 2 - ™V + " A)(m 2 - A 2 ) - A(P ) 2 )^} (4.50) 

The infinitesimal Lorentz boosts (4.49) can then be obtained by taking Poisson brackets 
with Ki defined in (4.28), 

«(A) = ^{m^ +BP »(A + (A2 " " A ' )} = {t.^'A',} , (4.51) 

in addition to (5p M = {p^, oj 0l Ki}. 

The evolution operator (4.16) can be applied to the relativistic system upon choosing Hq 
equal to p°. Then the mean values of x« and i evolve according to 

d A ~ 2 

^<*>> = -*<[^V = <|(l + |;)> 

^ < Li > = -i < [Up ] > + (^f ) = < (i + f,) > , (4.52) 



giving a mean velocity of 



d < Xi > < fH 1 + > 



(4.53) 



This is in contrast with the result in (3.25). So, as in the nonrelativistic case, the mean speed 
depends on the size and shape of the wavepacket as one approaches the noncommutative scale. 

Also as in the nonrelativistic case, the dynamical time operator (4.34) has continuous eigen- 
values which we again denote by r. Using the differential representation (2.4), the momentum- 
dependent eigenfunctions of t re \ are 

2 2 iA 2 (r-A) 

^A)=^^±^{^(m+pP)} m (4.54) 

The factor Aie lKT is determined from the orthonormality condition (4.26). After substituting 
in (4.54) and applying the integration measure (2.6), the left-hand side of (4.26) becomes 
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which after a change of integration variables from p to w = k + ^- ln|^-(m + p°)|, simplifies 
to 

r-co 

\M\ 2 / dw e i{T - T ' )w (4.56) 
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The delta function on the right-hand side of (4.26) is recovered for k — > -co and \M\ = 
1/(2V2ttA). 

4.3 ID conservative systems 

As an example of an interacting system we examine the simplest case of a one-dimensional 
conservative system. Following the procedure of the previous examples, we require the Hamil- 
tonian generating dynamics on Snyder space to be identical to the Hamiltonian Hq on canonical 
phase space, 

H = H = f - + V(x) , (4.57) 

and so the evolution parameter A cannot be identified with the time t. The relation between 
the two can once again be regarded as a gauge fixing condition. For this we should start with 
an action which is reparametrization invariant. Such an action can be written as 

<S C on = J d\[jj-V(x)i] , (4.58) 

with the dot again denoting differentiation in the evolution parameter A. It reduces to the stan- 
dard action for a conservative system when A equals t. The Hamiltonian constraint generating 
the reparametrization symmetry A — ^ A' (A) is now 

i = H + TT t w , (4.59) 

with p and 7Tj canonically conjugate to x and t, respectively. Since we are only considering one 
spatial dimension, the Snyder algebra in this case consists of the single relation between x and 
P, 

2 

{x,p} DB = l + ^ , (4.60) 

which we wish to have realized by the Dirac brackets. Then if *S>2 ~ denotes the gauge fixing 
condition its Poisson brackets should satisfy 

p 2 _ {a,$i}fep}-{g,Wi,p} -{#0,^2} f461) 

or simply 

2 

{H , ^2} = ^Tfj2 {*2, 7T t } (4.62) 

Explicit solutions for ^2 depend on the form for the potential energy V(x). For the example 
of a linear potential V(x) = —Fx, one has 

^ 2 = t+-^(Atan- 1 |-p) -A, (4.63) 
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up to an additive function of Hq. Setting ^2 (strongly) equal to zero, the evolution parameter A 
ranges over a finite domain of length irA/F for any classical trajectory, p = p so \ = Fi+ const ant. 
For such trajectories, A and t scale as in (3.14). The Dirac brackets of the time t with x and 
p resulting from the gauge fixing (4.63) are 

2 

{t, x}db = {*,p}db = , (4.64) 

which then can be promoted to commutation relations in the quantum theory. 

Another example is the one-dimensional harmonic oscillator with the standard Hamiltonian 
(3.8). Now (4.62) has the solution 

^2 = t + —{a tan" 1 amujx/p — tan" 1 mujxlp\ — A, a = — . = , (4.65) 

up to an additive function of Hq. This agrees with the result of the integration in (3.15) along 
classical trajectories x = x so \ = Asin(ujt + 4>), p = p so \ = mu)Acos(ujt + <j>) . Also, it reduces 
to the gauge fixing (4.7) for the free nonrelativistic particle in the limit oj — > 0. (This differs 
from the previous example where the free particle limit F — > was singular.) The resulting 
expressions for the Dirac brackets of the time t with x and p are 

{t,x}r>B = [mojxp — (A 2 + p 2 )atan _1 amux/p^ 

2 

{^p}db = ^ 4 (A 2 + p 2 ) yp + mojxa tan -1 amujx/pj (4.66) 

It is a nontrivial problem to find a quantization of this system. Also, the generalization of 
these examples to higher dimensions is not immediately obvious. 



5 Discussion 



In section 3, we gave conditions [Cf. (3.11)] for writing down any given dynamical system on 
the symplectic manifold associated with the Snyder algebra. Such a manifold is characterized 
by the symplectic two-form (3.1). Solutions, if they exist, may not be unique, and require fixing 
both the Hamiltonian H and the evolution parameter A in terms of the phase space variables 
and the time. Solutions were found for the case of free particle. There we could either set the 
time equal to the evolution parameter A, so that it remains a real parameter in the quantum 
theory, or in an alternative approach, equal to a monotonically increasing function of A and 
the phase space variables. In the latter, the time gets promoted to a Hermitean operator 
in the quantum theory. The standard energy-momentum dispersion relations were preserved 
in the latter approach, and the full Poincare (Galilei) symmetry group for the relativistic 
(nonrelativistic) particle was recovered. This was because the system could be derived from a 
Poincare (Galilei) invariant action. The introduction of a time operator in quantum theory has 
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nontrivial consequences, in particular with regards to causality violation. Although causality 
violation is a well-known feature of relativistic quantum mechanics, the situation here is more 
dramatic. From the evolution operator (4.16), we can write down a Green function which takes 
a state with time eigenvalue r', associated with evolution parameter A', to a state with time 
eigenvalue r, associated with evolution parameter A. It can be expressed as a function of r, r' 
and A — A' using the time eigenfunctions <fi T (p,\) [eqs. (4.25) or (4.54)], 

G(r, r', A - A') = J dfi(p) MP, A)*<M£ A') e'^'^ (5.67) 

The Green function does vanish when r — r' > and A — A' < 0, and so admits acausal 
time evolution in the quantum theory. Causality violation in relativistic particle mechanics is 
standardly cured by writing down the corresponding quantum field theory. However, this is 
unlikely to be the case here since the resulting noncommutative field theory is nonlocal. 

The question of whether or not there exist solutions to (3.11) when interactions are present 
can be nontrivial. In Section 3.2, we saw that for the case of one-dimensional systems we 
cannot in general set the evolution parameter equal to the time. So for general interactions, 
time is not associated with a commuting parameter in the quantum theory. In Section 4.3, 
we examined the specific case of an interaction with an external potential in one dimension. 
The action was written in a reparametrization invariant way and the time was obtained from 
a gauge fixing condition. Interactions with external potentials in higher dimensions should 
also be possible. Of particular interest are interactions between relativistic particles, which 
have been examined in a general context in [12]. A nontrivial co-product was introduced 
for the purpose of studying multiparticle states in a related work. [13] On the other hand, 
the treatment of particle interactions (in the absence of field theory) is problematic within 
our Hamiltonian framework, due to no interaction theorems. [14] Thus a reasonable direction 
for studying interactions for the systems discussed here is to develop field theory on Snyder 
space. Preliminary attempts to write down field theory on Snyder space were given in [15], [16], 
which rely on star product representations of the Snyder algebra. However, the proposed star 
products have certain inconsistencies, such as non-associativity. Since Snyder space is, in fact, 
a quantum lattice, it makes sense to instead consider a lattice field theory. Our work with 
particle dynamics on the quantum lattice, should serve as the ground work for writing fields 
on such a lattice. 
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